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Abstract 

The objective of this study is to verify the consistency of the pre- 
scription of alternative minimum coupling (connection) proposed by the 
Teleparallel Equivalent to General Relativity (TEGR) for the Dirac equa- 
tion. With this aim, we studied the problem of neutrino oscillations in 
Weitzenbock space-time in the Schwarzschild metric. In particular, we 
calculate the phase dynamics of neutrinos. The relation of spin of the 
neutrino with the space-time torsion is clarified through the determina- 
tion of the phase differences between spin eigenstates of the neutrinos. 
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1 Introduction 

The description of the gravitational field in the Teleparallel Equivalent to Gen- 
eral Relativity (TEGR) introduced by Maluf led to tensorial expressions for the 
energy, momentum and angular momentum of the gravitational field [T], [3J. 
This theory can be considered as a reformulation of Einstein's general relativity 
in terms of tetrad fields e a M , which is known as tetrad gravity. 

We examine the consistency of the Dirac equation in the TEGR through 
a new prescription of minimal coupling with the Dirac spinor fields ip. This 
alternative prescription, involves the Levi-Civita connection u) fl ab rather than 
spin connection of ui^ a b (field variable independent of the tetrad field e a M ). With 
respect to this connection, Maluf showed in 1994 [3] that it is possible to rule 
out LOjiab both in Lagrangian and Hamiltonian formulation. 

One motivation for this work, is the fact that, in the context of metric affinc 
theories of gravitation (MAG) , Obkuhov and Pereira [I] found an inconsistency 
in the coupling of the Dirac spinor gravitational field when the spin connection 
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w^afc was used in the covariant derivative. They found that the Dirac spinor 
fields couple with the gravitational field in a manner consistent only with spin 
or matter and with the spin tensor conserved. The TEGR was studied by 
Obkuhov and Pereira and can be considered to be a special case of the general 
theory of MAG. In 2003, Maluf showed that it is possible to bypass this prob- 
lem with the use of a new type of coupling to the Dirac spinor fields [S]. To 
check the coupling explicitly, we apply the Dirac equation of the TEGR with 
the connection °a> Ma & (totally dependent on tetrad field) to the problem of os- 
cillations of solar neutrinos or mixing in vacuum. In particular, we calculate 
the phase dynamics of neutrinos in space-time with torsion. The contributions 
of torsion and their relationship with the directions of spin in mass eigenstates 
were determined. These effects are compared with the structure of Minkowski 
and Riemann-Cartan geometry. In the Riemann geometry, there is no coupling 
between the spin of the particle and the gravitational field. 

The article is organized as follows: Section 2 presents a summary of the 
Hamiltonian formulation of the TEGR and the problem of using the spin con- 
nection ujfiab- In Section 3 we present the new minimum coupling to the Dirac 
equation in the TEGR and apply this equation to model neutrino oscillations 
in vacuum. The conclusions are presented in Section 4. 

Notation: Space-time indices /i,^, ... and local Lorentz SO(3, 1) indices 
a, run from to 3. Time and space indices are indicated according to fi = 0, 
i; a = (0), (i). The flat space-time is fixed by r] ab — e^e^g^ = (—+++). The 
tetrad field e a ^ and the arbitrary spin affine connection Lu^ ab yield the usual 
definitions of the torsion and curvature tensors: R a b^v = d^L0 ua b — d^tu^ab — 
T a — dfj,e a v+w^ - b e b v — ....The determinant of the tetrad field is represented 
by e = det(e a fj,). c is the speed of light, Ti is the Planck constant and G is the 
Newtonian gravitational constant. 



2 The Hamiltonian formulation of the TEGR 

In 1994, Maluf [3 implemented a Hamiltonian formulation of the TEGR with 
local symmetry by imposing Schwinger's time gauge 6J and found the following 
Hamiltonian density 



H(e u)i ,TlW\u kab ,P kab ,u; 0ab ) = -N k C k - NC + u 0a b J ab + (1) 

-Ne\ aMk R ablk - £ (m)(n) C7( m >^ + 
+d k (P kab uj ab) + ft [N k U kt + N(2eT i )] - 

where X ablk , N k , C (m *>^ J ab , A^and N are Lagrange multipliers, C k , S (ro)(n) , 
C, Ilr^i, J ab , Rabik are constraints given in Ref.[2|. P kab and EE"')* are compo- 
nents of the momentum canonically conjugate to io kab and eu)i, respectively. 
Maluf showed that the constraints in Eq. |T]) are first class by fixing uiQ ab = 0. 
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Then we can eliminate the momentum p kab from the Hamiltonian density. The 
evolution equation for LOkab leads to 



Ukab = {H,LO kab } =0, (2) 

and, for the equation to be consistent, it is necessary that tOkab — and then 
^fiat is ruled out of the theory. Then, the Lagrangian density and the field 
equations are invariant under global Lorentz transformations. 

3 Neutrino oscillations in the TEGR 

The tetrad field associated with the Schwarzschild metric was determined using 
Schwinger's time gauge [B] e^°'i = 0, e^) = and the condition for spatial 
symmetry g]: e (i)j = e (j)i 



( / \ 

o f- 1 

r 

\ rsinf? / 



( 3 ) 

\ rsinf? / 
where / = (1 — ^M-)^ , M is the mass of the sun. The determinant of the tetrad 



field is given by e = r 2 sin 6 



The equation proposed in the Maluf 's article [5] for the minimum coupling 
is described by a covariant derivative: 

i° 

= d^J - - w^E" 6 ^. (4) 



By substituting Jl} in the Dirac equation iTi^^D^ip — mocip = and using 
uab = "K^ ab where 
T v xn - Tpvx), we have: 



°<*Vl> = -Kfj, ab where K^ ab is the contortion tensor K Mab = \e a x e b v {T x ^ + 



itrf (9 M V + -K^^) - tnoc0 = 0. (5) 

By using the identity: e a bcf£ abcd = — 3!<5^,and by considering only the axial 
part of the contortion tensor represented by the vector A a — (A^ ,A), we 
obtain: 

K abc e abcd = e abcf Afe abcd = e abcf e abcd A^ = -3!#AA (6) 
By using the identities involving the Dirac matrices 7° = e a M 7 M = , 7W): 

J a ll b , 71 = -2v a "j c + 2v ac J b + 2ie afccrf 7( 5 )7 (0) 7 (1) 7 (2) 7 (3) 7d, ( 7 ) 
7(5) = _ i7 (0) 7 (D 7 (2) 7 (3) 5 (8) 

£ ab = ^[7 a ,7 fc ], (9) 
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7 (0) 7(5)7a 



a A a = 1(5) A^ + t-l, 



7( 

if 
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7 Q 7 b + 7 V 



and by denning the components of the momentum as 
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we can write the Dirac Hamiltonian matrix with the help of Eq. [6] 



(10) 

(11) 
(12) 

(13) 
(14) 
(15) 



H = 



( fIm c 2 + pcf(a^A^ + 

f 2 Ia^p r c + fIaWp gC + 
V +fIa^p„c-pcfA^I 



f 2 IaWp r c+fIaWp e c+ 
+fI<rWp v c-%hcfAWl 
-flm c 2 + pcf(a^AW + 
+a^A^+a^A^) 



J 



(16) 

where <jw are the Pauli matrices, and I is the identity matrix. 

First, we analyze the azimuthal motion which, although is not realistic, is 
more simple and was considered in Ref. [7] and [5] in a different formalism. 

3.1 Azimuthal Motion 

Zhang and Pereira in Ref. [9] was one of the first to calculate the mass neutrino 
oscillation induced by torsion in the context of the weak gravitational expansion 
method. For the calculation of the phase difference in azimuthal and radial 
motion, we proceed in the same way as in Ref. [10] . The Dirac Hamiltonian for 
azimuthal motion is given by (|16p . where we have substituted ~]f — (p r ,pe,p v ) — 
(0,0, p). Then the positive eigenvalue of Eq. (Til)]) i s: 



E+=p{l 



2cV 



3ha^A^\ . 3 fc _ 

)cf + -he fa ■ A, 

2p ' J 2 J 



where we used the ultra-relativistic limit E ~ cp 3> mc 2 . The negative engen 
value gives similar results. 



With the help of the unitary transformation for the spinor ij) = U 
where U is a the unitary matrix, we can write 
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E+Z+= [fpc+ f -^-^U + + (17) 

+ \ lhcf(A^+iAW) -%hcf (A® - A&) )Z + - tn W+> 

where we found the following eigenvalues of Eq. p7|) for the mass eigenstates 
with spin up and spin down 



/jK+ M£! + | Bo/A | s ; = iR » £ t, (18) 

respectively. Here A = y/(AM - A< 3 )) 2 + (At 1 )) 2 + (A^) 2 . 

The phase of the neutrino is obtained by the integration $ = i J E+dt, 
resulting in 

♦ I ( fER A fm 2 c 3 R^ 3„,„„ A 1 

h\c (f+ IE A( P+2 nfARA Vf- ( 19 ) 



*i = - [ ^ + J -^-*p ~ - 2 nfARA<p\ 



2„3 j 

(20) 



for spin up and spin down, respectively. We used the ultra-relativistic limit 
pc ~ E and cdi ~ Rdtp, where R is the radius for the circular orbit of the 
neutrino and Aip is the angular dislocation. 
Now we have three possibilities 

1. From Eqs. (p~9j) and ([20]) . the phase difference between eigenstates that 
have the same spin state is given by 

(21) 

2. For the first mass eigenstate with spin down and the second eigenstate 
with spin up, we obtain 

^=^-^={mm-' fA } RA ^ (22) 

3. For the first mass eigenstate with spin up and the second eigenstate with 
spin down, we obtain 

A ^^-^={mrn +3fA } RA ^ (23) 

where Am 2 = nvk — m? . 
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3.2 The radial motion 

For the calculation of the phase difference in radial motion, we proceed in the 
same way as in the case of azimuthal motion. Substituting ~jt = (p r ,Pe,Ptp) = 
(p, 0, 0) in the expression (|f 6[) , and using the ultra- relativistic aproximation pc ~ 
E, we obtain the following phase differences for the neutrino mass eigenstates: 

1. For eigenstates with the same spin: 

A$ = $| - $J = ^H-^Ar. (24) 
2 1 2Eh V ' 

where Ar = rs — ta- Here r a is the radius of the sun where the neutrinos was 
produced and tb is the distance from the center of the sun to the earth's surface 
where the neutrino was measured. 

2. For the first eigenstate with spin down and the second eigenstate with 
spin up 

A $ = - $T = ^ff A r - 3A(Ar - M° In % (25) 

Zhjfl c T A 

3. For the first eigenstate with spin up and the second eigenstate with spin 
down 

A$ = $^-$| = — Ar + 3A(Ar —In—), (26) 

IhjTi c tj± 

where A = y/(AM - A^) 2 + {A^) 2 + (A^) 2 . Here, we used cdt ~ dr, f = 

\ rc^ ' rc z 

Note that the magnitude of the contribution from the torsion is very small 
compared with the mass differences between neutrinos [TU] . 

If electron neutrinos are produced at time t = 0, the time evolution of the 
flavor eigenstate v e is given by: 

v e (t) = cosee-** 1 ^! +sinee^* 2(t) z/ 2 , (27) 

where v\ and v-i ae called mass eigenstates at t — and 8 is called the "mixture" 
angle. The probability of measuring a muon or tau ineutrino is given by: 

P{v e -> v^ T ) = sin 2 26 sin 2 (28) 



4 Conclusions 

In azimuthal and radial motion, the torsion makes a contribution to the phase 
dynamics of neutrinos that depends on the directions of spin in the mass eigen- 
states. If both mass eigenstates have the same direction of spin, there is no 
contribution to oscillation from torsion, but if the mass eigenstates have oppo- 
site spin directions there is a contribution to the oscillation of neutrinos that 
comes from torsion.. 
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Without torsion or curvature (with / = 1), our results coincide with results 
for neutrino oscillations in vacuum in flat space-time [7]. 

Our result formally coincides with a result obtained by Adak et al. [TU], from 
the Dirac equation in a Riemann-Cartan geometry with the connection of spin 
k-Vaf> — k-Vafc + K/j,ab- However, in radial motion, the radial component of mo- 
mentum p r found in our paper is less than the radial component of momentum 
found by Adak et al. Therefore, we can claim that the energy of a neutrino in 
the space-time of a teleparallel geometry, called Weitzenbock space-time, is less 
than the energy of a neutrino in Riemann-Cartan space-time and the probability 
of transition to a muon or tau neutrino is higher in the teleparallel geometry. 

The prescription of minimum coupling adopted by Maluf does not lead to 
inconsistencies of the type that Obukhov and Pereira found and are qualitatively 
consistent with the literature, at least for the problem of neutrino oscillation. 
The formal inconsistency of the Dirac equation is removed if we adopt the Levi- 
Civita connection of the theory as uj tia b = —K^ a b leaving the total Lagrangian 
density in the TEGR in the presence of Dirac spinor fields invariant under global 
Lorentz transformations. 
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